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Claims : 

1 1. A method of processing a digital image, comprising the 

2 steps of: 

3 providing digital data indexed to represent positions of 

4 an image having S spectral bands for simultaneous output on a 

5 display, said digital data being indicative of an intensity 

6 value I-(x,y) for each position (x,y) in each i-th spectral 

7 band; 

8 defining a classification of said image based on dynamic 

9 range of said image in each of said S spectral bands; 

10 adjusting said intensity value for said each position in 

5| each i-th spectral band to generate an adjusted intensity 

S| value for said each position in each i-th spectral band in 

fft- 

f3 accordance with 

^W D (log I d (x,y) - log[I i (x t y) *F n (x,y) ] ) , 2 = 1,..., S 

111 where S is the number of unique spectral bands included in 

£3; said digital data and, for each n, W n is a weighting factor 

16 and F n (x,y) is a unique surround function applied to said each 

17 position (x,y) and N is the total number of unique surround 

18 functions; and 

19 filtering said adjusted intensity value for said each 

20 position of said image in each of said S spectral bands using 

21 a filter function based on said classification of said image 

22 wherein a filtered intensity value Rj(x,y) is defined. 
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2. A method according to claim 1 wherein each said unique 
surround function is a Gaussian function . 
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3. A method according to claim 2 wherein said Gaussian 
function is of the form 



e * 



satisfying the relationship 



^4 



k n j J e c " dxdy = 1 



where 



S5 



L*6 



r = \/x 2 +y 2 

and, for each n, k n is a normalization constant and c n is a 
unique constant for each of said N unique surround functions. 
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4. A method according to claim 1 further comprising the step 
of multiplying said filtered intensity value R-(x,y) by 



log 



Bl^x.y) 
s 

i 

i=i 



to define a color-restored intensity value R'^x^y), where B 
is a constant - 
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5. A method according to claim 1 wherein said each position 
(x,y) defines a pixel of said display. 
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1 6. A method according to claim 1 wherein, for each n, W n =l/N. 

1 7. A method according to claim 1 wherein said step of 

2 defining comprises the step of using image statistics 

3 associated with said image in each of said S spectral bands to 

4 select said filter function. 

1 8. A method according to claim 7 wherein said image 

2 statistics include brightness and contrast of said image in 

3 each of said S spectral bands. 

-3 ; 

y| 9. A method according to claim 1 further comprising the steps 

y[ selecting a maximum intensity value V^^y) from the 

group consisting of said intensity value I-(x,y) and said 

jjf filtered intensity value R f (x,y); and 

Ml displaying an improved image using said maximum intensity 

H value V i (x, y) . 

1 10. A method according to. claim 4 further comprising the 

2 steps of: 

3 selecting a maximum intensity value V^x^y) from the 

4 group consisting of said intensity value I 1 - (x ; y) and said 

5 color-restored intensity value R'^Xjy); and 

6 displaying an improved image using said maximum intensity 

7 value V 1 - (x, y) . 
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11. A method of processing a digital image, comprising the 
steps of: 

providing digital data indexed to represent the positions 
of a plurality of pixels of a J-row by K-column display, said 
digital data being indicative of an intensity value I(x,y) for 
each of said plurality of pixels where x is an index of a 
position in the J-th row of said display and y is an index of 
a position in the K-th column of said display wherein a JxK 
image is defined; 

convolving said digital data associated with each of said 
plurality of pixels with a function 

-r 2 

e c 

to form a discrete convolution value for each of said 
plurality of pixels, said function satisfying the relationship 



k is a normalization constant and c is a constant; 

converting, for each of said plurality of pixels, said 
discrete convolution value into the logarithm domain; 

converting, for each of said plurality of pixels, said 
intensity value into the logarithm domain; 
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20 subtracting, for each of said plurality of pixels, said 

21 discrete convolution value so-converted into the logarithm 
2 2 domain from said intensity value so-converted into the 

23 logarithm domain, wherein an adjusted intensity value is 

24 generated for each of said plurality of pixels; and 

25 filtering said adjusted intensity value for each of said 

26 plurality of pixels with a filter function that is based on 

27 dynamic range of said JxK image wherein a filtered intensity 

28 value R(x,y) is defined. 

Jf: 12. A method according to claim 11 wherein the value of said 

Cg: constant c is selected to be in the range of approximately 

■fa. is. . , 

%l 0.01 to approximately 0.5 of the larger of J and K. 

f\ 13. A method according to claim 11 further comprising the 
steps of : 

Jff selecting, for each of said plurality of pixels, a 

maximum intensity value V(x,y) from the group consisting of 

5 said intensity value I(x,y) and said filtered intensity value 

6 R (x, y ) ; and 

7 displaying an improved image using said maximum intensity 

8 value V(x,y) . 
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1 14. A method of processing a digital image, comprising the 

2 steps of: 

3 providing digital data indexed to represent the positions 

4 of a plurality of pixels of an J-row by K-column display, said 

5 digital data being indicative of an intensity value I f (x,y) 

6 for each i-th spectral band of S spectral bands for each of 

7 said plurality of pixels where x is an index of a position in 

8 the J-th row of said display and y is an index of a position 

9 in the K-th column of said display wherein a ( JxK) 1 image is 
10 defined for each of said S spectral bands and a JxK image is 
14 defined across all of said S spectral bands; 

M defining a classification of said JxK image based on 

ffl dynamic range of each said (JxK) • ; 

m 

1=4 convolving said digital data associated with each of said 

i§ plurality of pixels in each i-th spectral band with a function 



e c " 



j| for n=2 to N to form N convolution values for each of said 

17 plurality of pixels in each said i-th spectral band, said 

18 function satisfying the relationship 



k n J J e c " dxdy = 1 



19 where 



i = \/x 2 +y 2 

20 and, for each n, k n is a normalization constant and c n is a 
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21 unique constant; 

22 converting, for each of said plurality of pixels in each 

23 said i-th spectral band, each of said N convolution values 
2 4 into the logarithm domain; 

25 converting, for each of said plurality of pixels in each 

26 said i-th spectral band, said intensity value into the 
2 7 logarithm domain; 

28 subtracting, for each of said plurality of pixels in each 

29 said i-th spectral band, each of said N convolution values so- 
il converted into the logarithm domain from said intensity value 

4% 

1J so-converted into the logarithm domain, wherein an adjusted 

nj 

If intensity value is generated for each of said plurality of 

»| pixels in each said i-th spectral band based on each of said 

34 N convolution values; 

3 

^| forming a weighted sum for each of said plurality of 

3j€ pixels in each said i-th spectral band using said adjusted 

© intensity values; and 

38 filtering said weighted sum for each of said plurality of 

39 pixels in each said i-th spectral band with a filter function 

40 that is based on said classification of said JxK image wherein 

41 a filtered intensity value R^X/y) is defined. 

1 15. A method according to claim 14 wherein the value for each 

2 said unique constant c n is selected to be in the range of 

3 approximately 0,01 to approximately 0.5 of the larger of J and 

4 K. 
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16. A method according to claim 14 further comprising the 
step of multiplying said filtered intensity value R^x^y) by 

[ Bl^x.y) 
log 

. i = l 

to define a color-restored intensity value R'^x^y), where B 
is a constant and S is a whole number greater than or equal to 
2. 

17. A method according to claim 14 wherein said step of 
defining comprises the step of using image statistics 
associated with each said ( JxK) i image to select said filter 
function . 



18. A method according to claim 17 wherein said image 
statistics include brightness and contrast of each said (JxK) • 
image . 




1 19. A method according to claim 14 further comprising the 

2 steps of: 

3 selecting a maximum intensity value (x,y) from the 

4 group consisting of said intensity value 1^ (x,y) and said 

5 filtered intensity value R-(x,y); and 

6 displaying an improved image using said maximum intensity 

7 value V i (x, y) . 

1 2 0. A method according to claim 16 further comprising the 

^2 steps of: 

y3* selecting a maximum intensity value V-(x,y) from the 

55 group consisting of said intensity value I^x^) and said 

\§ color-restored intensity value R'jfx^) ; and 

yS displaying an improved image using said maximum intensity 

H7 value V- (x, y) . 



m 
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1 21. A method of processing a digital image, comprising the 

2 steps of: 

3 providing digital data indexed to represent positions of 

4 an image having S spectral bands for simultaneous output on a 

5 display, said digital data being indicative of an intensity 

6 value I^Xjy) for each position (x,y) in each i-th spectral 

7 band; 

8 defining a classification of said image based on dynamic 

9 range of said image in each of said S spectral bands; 

10 adjusting said intensity value for said each position in 

IJ each i-th spectral band to generate an adjusted intensity 

|3> value for said each position in each i-th spectral band in 

IB accordance with 

B £ W n (log I ± {x,y) - loqll^x.y) *F n (x,y) ) ) , i = l,...,S 

|2 where S is a whole number greater than or equal to 2 and 

MS defines the total number of spectral bands included in said 

•SKIS.. , • _ 

it digital data and, for each n, W is a weighting factor and 

17 F n (x,y) is a unique surround function of the form 

-r 2 

e ^ 

18 satisfying the relationship 

-r 2 

k n f fe c?n dxdy = 1 

19 where 
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r = >Jx 2 +y 2 

20 and, for each n, k n is a normalization constant and c n is a 

21 unique constant where N is the total number of unique surround 

22 functions ; 

23 filtering said adjusted intensity value for said each 

24 position in each i-th spectral band with a function based on 

25 said classification of said image wherein a filtered intensity 

26 value R^x^y) is defined; and 

27 multiplying said filtered intensity value R f (x,y) by 
O log 

fl to define a color-restored intensity value R'^x^y) , where B 

21) is a constant. 

3 

gi 22. A method according to claim 21 wherein, for each n, 

ft i 

III W n =l/N. 

1 23. A method according to claim 21 wherein the value for each 

2 said unique constant c n is selected to be in the range of 

3 approximately 0.01 to approximately 0.5 of the larger of J and 

4 K. 
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1 24. A method according to claim 21 wherein said step of 

2 defining comprises the step of using image statistics 

3 associated with said image in each of said S spectral bands to 

4 select said filter function. 

1 25. A method according to claim 24 wherein said image 

2 statistics include brightness and contrast of said image in 

3 each of said S spectral bands. 

f? 4 26. A method according to claim 21 further comprising the 

?S steps of: 

selecting a maximum intensity value V,- (x,y) from the 

Nji group consisting of said intensity value I.(x,y) and said 

M5 color-restored intensity value R'^x^) ; and 

y5 displaying an improved image using said maximum intensity 

f|J7 value V i (x, y) . 
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